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Abstract. The basic notion for a motion of a heavy rigid body fixed 
at a point in three-dimensional space as well as its higher-dimensional 
generalizations are presented. On a basis of the Lax representation, 
the algebro-geometric integration procedure for one of the classical 
cases of motion of three-dimensional rigid body - the Hess-Appel'rot 
system is given. The classical integration in Hess coordinates is pre- 
sented also. For higher-dimensional generalizations, the special at- 
tention is paid in dimension four. The L-A pairs and the classical in- 
tegration procedures for completely integrable four-dimensional rigid 
body so called the Lagrange bitop as well as for four-dimensional 
generalization of Hess-Appel'rot system are given. An n-dimensional 
generalization of the Hess-Appel'rot system is also presented and its 
Lax representation is given. Starting from another Lax representa- 
tion for the Hess-Appel'rot system, a family of dynamical systems 
on e(3) is constructed. For five cases from the family, the classi- 
cal and algebro-geometric integration procedures are presented. The 
four-dimensional generalizations for the Kirchhoff and the Chaplygin 
cases of motion of rigid body in ideal fluid are defined. The results 
presented in the paper are part of results obtained in the last decade. 
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1. Introduction 

The rigid body motion is one of the most studied and most interesting systems 
of classical mechanics. Nevertheless, there are still some important open questions 
and problems concerning it. In this paper we will focus on the problem of integra- 
bility of motion of heavy rigid body fixed at a point. We are going to present the 
classical integration of some known integrable cases as well as the algebro-geometric 
integrations procedure based on existence of the Lax representation. 
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One of the main questions in a study of the system of differential equations of 
motion of some mechanical system is integrability or solvability. The notion of 
intcgrability is very close to the existence of the first integrals, i.e. functions that 
are constants on solutions of the system. Until the beginning of XX century, the 
theory of integrable system had been intensively developed with great influence of 
leading mathematicians and mechanicians of that time (Euler, Hamilton, Jacobi, 
Lagrange, Poincare, Liouville, Noether, Kowalevski an many others). For proving 
integrability they usually used some of the basic methods of that time: method of 
separation of variables and Noether's theorem for finding integrals from symmetries. 
It became clear that algebraic geometry and theory of theta functions, that was 
intensively developed in that time, have an important role in integration of the 
dynamical systems. For example, the solutions of the Euler and Lagrange cases of 
motion of a rigid body fixed at a point are meromorphic functions on an elliptic 
curve. Starting from that fact, Sofia Kowalevski formulated the problem of finding 
all cases of rigid body motion fixed at a point whose solutions are unique functions 
of complex time that admit only moving poles as singularities. She proved that this 
is possible only in one more case, today called the Kowalevski case. She found the 
additional first integral and she completely solved the system in theta functions. 
The importance of the Kowalevski paper is reflected in the fact that thousands of 
papers are devoted to it. For recent progress, geometric interpretation and certain 
generalizations of the Kowalevski top see jT5J [21] • 

In the 60's of XX century the big progress was made in the theory of inte- 
grable systems. It was proved that some nonlinear partial differential equations 
(Korteveg-de Vries (KdV), Kadomtsev-Petviashvili (KP) and others) are infinitely- 
dimensional Hamiltonian systems. Also, a new method appeared: algebro-geometric 
integration procedure. It is based on the existence of a Lax representation (or L-A 
pair). A system admits L-A pair with spectral parameter if there exist matrices 
L(X),A(X) such that equations of the system can be written in the form: 

(1.1) j t L(X)=[L(X),A(X)], 

where A is a complex number. An important case, when L(X) and A(X) are matrix 
polynomials in A, was studied by Dubrovin in [25] (see also HZ])- The first 
consequence of (|1.1[) is that the spectrum of matrix L(X) is a constant function in 
time, i.e. coefficients in spectral polynomial are first integrals. If, from L-A pair, 
one gets enough integrals for integrability, then the system can be integrated using 
algebro-geometric integration procedure, which is developed by the Novikov school. 
In that procedure, the Baker-Akhiezer function plays the key role. This function is 
common eigenfunction of operators 4r+A(X) and L(X), defined on the spectral curve 
r naturally associated to L-A pair. The Baker-Akhiezer function is meromorphic 
on T except in several isolated points where it has essential singularities. For a 
detailed explanation see [25j [2H [23 [HI E [M] ■ Let us mention also that the Lax 
representation is useful for constructing higher-dimensional generalizations of the 
system. In pQ Adler and van Moerbeke have presented an additional approach 
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for integrability. Both methods have been successfully applied to the rigid body 
motion (see 02 Q3J H E] ) . 

The theories of rigid body motion and of intcgrable dynamical systems have 
been intensively studied by Serbian scientists (see for example books and mono- 
graphs [TUl H3 HI E] ) • At the Seminar Mathematical Methods of Mechanics in the 
Mathematical Institute SANU, supervised by Vladimir Dragovic, a group of young 
researchers including myself, Milena Radnovic and Bozidar Jovanovic, the theory 
of integrable dynamical systems has been studied for almost 20 years. Here we will 
review some of the joint results obtained with Vladimir Dragovic in the last decade 

(see [ED USES EH |22]). 

This paper is organized as follows. In Section 2 the notions of Poisson struc- 
ture and integrability in Liouville sense are given. Also the basic steps in algebro- 
geometric integration procedure are performed. The basic facts about three-dimensional 
motion of a rigid body are presented in Section 3. The classical as well as the 
algebro-geometric integration procedures for the Hess-Appel'rot case of motion of 
three-dimensional rigid body are given also. The basic facts on higher-dimensional 
rigid body motion as well as the definition of the Lagrange bitop and n-dimensional 
Hess-Appel'rot systems are presented in Section 4. In Section 5 we present a con- 
struction of a class of systems on the Lie algebra e(3). For the five cases when 
an invariant measure is preserved, the classical and algebro-geometric integration 
procedures are given. The four-dimensional generalizations of the Kirchhoff and 
Chaplygin cases of the motion of the rigid body in ideal fluid are given in Section 
6. 

2. Poisson structure on manifolds. Integrability. Algebro-geometric 

integration procedure 

Let M be a smooth manifold, and C°°(M) algebra of smooth functions on M. 

Definition 2.1. A Poisson bracket on M is a map { , } : C°° (Af) x C°° (M) -> 
C°°{M) that for f,g,he C°°{M) satisfies: 

(1) bilinearity: {A/ + fig, h} = \{f, h} + fi{g, h}, A,/ieM 

(2) skew-symmetry: {f,g} = -{g,f} 

(3) Leibnitz rule: {fg, h} = f{g, h} + g{f, h} 

(4) Jacobi identity: {/, {g, h}} + {g, {h, /}} + {h, {/, g}} = 0. 

If (x ,x , ...,x n ) are coordinates on M, using the Leibnitz rule one has: 

where P %3 = {x l , x 3 }. Poisson bracket is also called Poisson structure on manifold, 
and a manifold endowed with Poisson structure is a Poisson manifold. 

Poisson bracket can be degenerate. Then matrix P = (P* J ) is singular. If P 
is nonsingular, then, because of skew-symmetry, the dimension of M is even and 
inverse matrix P~ l gives symplectic structure on M. Functions whose Poisson 
bracket with any other function is equal to zero are called the Casimir functions. 
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For a smooth function H on manifold M, the system of equations 

= = If 

j 

is called Hamiltonian system with the Hamiltonian function H. The vector field 
Xjj = Y^j P lJ §§r is called the Hamiltonian vector field associated with H. 
A function / is a first integral of a system of differential equations 

x l = X l (x), i = 1, n 

if it is constant along every solution of the system, or in other words, if J2i §^^ 1 — 
X{f) = 0. Geometrically, it means that each solution lies on a hypersurface / = 
const. For the integrability in quadratures one needs n — 1 first integrals. However, 
by the Jacobi theorem, if a system preserves the standard measure, i.e. if the 
divergence of the vector field X is zero, then for the integrability in quadratures 
one needs only n — 2 first integrals (see [Jj 133] ). 

For Hamiltonian systems there is an additional structure, the Poisson structure. 
A function / is a first integral of a Hamiltonian system with the Hamiltonian 
function H if it Poisson-commutes with the Hamiltonian, i.e. if {/, H} = 0. The 
following theorem is fundamental concerning the integrability of the Hamiltonian 
systems. 

Theorem 2.1 (Liouville-Arnol'd). Let M 2n be a symplectic manifold, and f\ — 
H, /2, f n functions that satisfy {fi,fj} = 0. Denote Mf — {x £ M\fx(x) = 
const, f n (x) — const.}. If fx, f n are functionally independent on Mf then 

(1) Mf is a smooth manifold invariant under the Hamiltonian flow with the 
Hamiltonian H = f\. 

(2) // Mf is compact and connected, then it is diffeomorphic to n-dimensional 
torus T n . 

(3) There exist coordinates if — {ifi, ip n ) on T n in which the Hamiltonian 
flow is linearized: 

<Pi =Wj(/i, ...,/„). 

(4) The Hamiltonian system with Hamiltonian H can be solved in quadratures. 
For proof see [B] . 

The Hamiltonian system that satisfies the Arnol'd-Liouville theorem will be 
called completely integrable in the Liouville sense. 

2.1. The basic steps of the algebro-geometric integration procedure. We will 

give here a short description of algebro-geometric integration procedure. For details 

see [23 [IE [271 [23 [El 0. 

The existence of a Lax representation for a system of differential equations 
is equivalent to commutativity of operators: 



j t +A(t,X),L(t,X) 



0. 
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Let <E>(i, A) is the fundamental solution matrix for the equation: 

(2.1) (| + A(t,A))*(t,A)=0, 

normalized by the condition $(0, A) = 1. From the Lax representation one gets 
that L(t, A)$(i, A) is also a solution of (|2.1[) . Since every solution is determined by 
its initial conditions, we have 

L(t,A)$(i,A) = $(t,X)L(0,X). 

Consequently, the matrices L(t, A) and L(0, A) have the same spectrum. In other 
words, the coefficients of the characteristic equation 

(2.2) p(X, n) = det(L(A) - \i ■ 1) = 

are the first integrals of the system. 

The equation (I2.2[) defines algebraic curve F, called the spectral curve. Let L 
and A are n X n matrices. Over A € C we have n points on T with coordinates 
(A,/ii), (A,/i n ). To each of these points corresponds eigenvector h(t,(X,)i)) of 
matrix L(t, A): 

L{t, \)h(t, (A, jUfc)) = n k h(t, (A, Hk)). 

Fix the normalization 

n 

Y / h l {t,(X^)) = l. 

i=l 

Normalized vector h can be regarded as meromorphic vector- function on T. Intro- 
duce the function 

^(t,(A,/z)) = $(i,A)/i(0,(A,/i)). 

In what follows we will see that this function has a key role in the integration 
procedure. 

On can easily check that ip(t, (X, /j,)) satisfies the following relations: 
L(t, X)ip(t, (A, /i fe )) = /ifeV(*, (A, jtxfc)), 

(^+A(t,A))^,(A, Mfc )) = 0. 

The following theorem gives us the analytical properties of ip(t, (A, //)). Denote 
with Pi, P n the points over A = oo. 

Theorem 2.2. [27] 77ie vector-function ip(t, (A, /x)) /ias f/ie following properties: 

(1) It is meromophic on T\{P\, P n }. Its divisor of poles has degree g+n—1, 
and it does not depend on time, where g is the genus of the curve T. 

(2) In the neighborhood of Pk the function ip(t, (A, //)) has the form: 

i>{t, (A, A*)) = «(Mfc ^exp^t, z k )}, 

where z^ 1 is a local coordinate in the neighborhood of Pk, a is the holo- 
morphic vector-function, and qk are polynomials. 
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The functions that satisfy conditions from thcorcni l2.2l arc called n-point Baker- 
Akhiezer functions. The example of such function is the exponential function e 1 / 2 
on CP 1 . It is holomoprphic everywhere except in point z = 0, where it has an 
essential singularity. The most general definition of the Baker- Akhiezer functions 
is given by Krichever. For a history and details see [27] . 

Definition 2.2. (see [27]) Let Pi,...,P n are points on a Riemann surface T of 
genus g and z^ 1 are local coordinates in the neighborhoods of these points such 
that z^ 1 {P k ) = 0. Let qi(z), qk(z) are polynomials and let V is divisor on F. 
The n-point Bakcr-Akhiezer function ip is a function that satisfies the following 
conditions: it is meromorphic on F \ {P%, P n }, its divisor of zeros and poles 
satisfies (ip) + T> ^ 0, and in the neighborhood of the each point Pk the function 
ip(P) exp(-q k (z k (P))) is analytic. 

The basic idea of algebro-geometric integration procedure is to reconstruct ij;(t, (A, /z)). 
From the analytical properties given in theorem l2.2l one can in terms of theta func- 
tions explicitly construct ip(t, (A,/z)), and using it solve the system. 

Theorem 2.3. For a non-special divisor T> of degree N the dimension of the linear 
space of functions with properties from definition \2.2\ is equal N—g+1. Particularly, 
if the degree of T> is g, then ip is uniquely determined up to factor and it is given 
by: 

^hi' p ' 6(A{T) - A(V) - K) 

where VL qk are normalized Abelian differentials of the second order, which in the 
neighborhood of Pk has the form 

£l qk = dqk(zk(P)) + holomorphic part, 

the vector j s a vector of b-periods of the differential Q qk , and K is the vector 

of Riemann constants. 

As a corollary, in a generic situation, one has n functions vj) 1 , ...,ip n that sat- 
isfy the theorem 12.21 Let ^(t, A) be the matrix of which the columns are vectors 
tjj(t, (A, Ufe)),fc = 1, ...,n (here /Xfc are eigenvalues of £(A)). One has 

L(t,\) = *fM- 1 , A(X) = ^-\ 

where fi = diag(/ii, Consequently, in a general situation, from the Baker- 

Akhiczcr functions one can find the matrix L and A as functions of time, or in 
other words one can integrate the system. 

3. Motion of a heavy rigid body fixed at a point 

A three-dimensional rigid body is a system of material points in R 3 such that 
the distance between each two points is a constant function of time. We shall 
consider motion of a rigid body with a fixed point O. The configuration manifold 
is the Lie group 5*0(3). Two different Euclidian coordinate frames are associated 
to the system: the first one Oxyz is fixed in the space, and the second, moving, 
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OXYZ is fixed in the body. With the capital letters we will denote elements of the 
moving reference frame, while the lowercase letters will denote elements of the fixed 
reference frame. If some point of the body has the radius vector Q in the moving 
coordinate system, then its radius vector in the fixed frame is q(t) = B(t)Q, where 
B(t) <E SO(3) is an orthogonal matrix. The velocity of that point in the fixed 
reference frame is given by 

v(t) = = B(t)Q = B{t)B-\t)q{t) = w(t)flt), 

where uj{t) = BB~ X . It can be proved that uj is an skew-symmetric matrix. Using 
the isomorphism of (]R 3 , x), where x is the usual vector product, and (so(3), [, ]), 
given by 



(3.1) a = (ai, a 2 , a 3 ) i-> a 



-a 3 a 2 
03 — ai 
— a,2 &i 



matrix u>(t) is corresponded to vector uj(t) - angular velocity of the body in the 
fixed reference frame. Then v(t) — uj(t) x q{t). One can easily see that uj(t) is the 
eigenvector of matrix uj(t) that corresponds to the zero eigenvalue. 

In the moving reference frame, V(t) = B(t)~ 1; v(t), so V(t) — f2(t) x Q, where 
£l(t) is the angular velocity in the moving reference frame and corresponds to the 
skew-symmetric matrix f2(t) = B~ 1 {t)B{t). 

At the fixed moment of time, the vector Q defines the line I through the fixed 
point. For points on I, vectors q(t) and uj(t) are collinear. Hence, velocities of 
points on I are equal to zero. The line I is called instantaneous axis of rotation. 

The existence of instantaneous axis of rotation can be regarded as a infinitesimal 
version of Euler's rotation theorem. The theorem states that any finite displacement 
of a rigid body with fixed point is equivalent to the rotation about some axes 
through the fixed point. Here we will present the original proof given by Euler in 
1775 J25]. 

Theorem 3.1 (Euler 29 ). If a sphere is turned about its centre, it is always 
possible to assign a diameter, whose direction in the displaced position is the same 
as in the initial position. 

Proof. Consider a great circle k± in initial state, which after the displacement goes 
to a circle k 2 . Denote by A a point of the intersection of these two circles. Since 
A belongs to k\, after rotation it goes to the point a g k 2 . On the other hand, the 
point A is on k%, so there is point a G k% which after rotation goes to A. We will 
show that there exists a point 0\ on the sphere, which equally refers to the circle 
fci as to the circle k 2 . If we suppose that 0\ is constructed, then the arcs 0\A 
and 0\a should be equal to each other. Also, the arcs 0\A and 0\a are similarly 
inclined towards the circles k\,k 2 . Consequently, the angles 0\aA and 0\Aa are 
equal also. But since the arcs 0\a and 0\A are equal, the angles 0\aA and 0\Aa 
are also equal, whence 0\Aa — 0\Act. It is clear that 0\ lies on the arc bisecting 
the angle aAa. So, 0\ can be constructed as an intersecting point of the sphere 
and the following two planes through the fixed point O (the center of the sphere) . 
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The first plane is the symmetry plane of the angle aAa and the second one is the 
symmetry plane of the arc Aa. The diameter that we are looking for is determined 
by the point Ox and the center of the sphere. □ 




Let us stress that it is natural to consider the angular velocity as an skew- 
symmetric matrix. The element W12 corresponds to the rotation in the plane de- 
termined by the first two axes Ox and Oy, and similarly for the other elements. 
In the three-dimensional case we have a natural correspondence given above, and 
one can consider the angular velocity as a vector. But, in higher-dimensional cases, 
generally speaking, such a correspondence does not exist. We will see later how in 
dimension four, using isomorphism between so(4) and so(3) x so(3) two vectors in 
the three-dimensional space are joined to an 4 x 4 skew-symmetric matrix. Since 
we cannot imagine higher-dimensional world (or at least, it is not easy to imagine 
it), it is much easier to consider, for example, two dimensional world. If the two- 
dimensional people consider rotation of a rigid body with fixed point, they conclude 
that angular velocity is a two-dimensional skew-symmetric matrix. If the element 
UJ12 is positive, then rotation goes in the positive sense, otherwise it goes in the 
negative one. There is no third dimension, so they cannot conclude that it can be 
seen as a rotation about z-axis. In the two-dimensional world the z-axis does not 
exist! 



10 



BORISLAV GAJIC 



The moment of inertia with respect to the axis u, defined with the unit vector 
u through a fixed point O is : 



where di is the distance between z-th point and axis u, and I is inertia operator 
with respect to the point O defined with 



In coordinates (X,Y,Z) the diagonal elements h\,l22,hz 01 I arc moments of 
inertia of the body with respect to the coordinate axes OX, OY, OZ respectively. 
For example In — X) m i(^i 2 + Non-diagonal elements are called centrifugal 
moments of inertia. For example I12 = — J] fn%^-iXi and similar for other Jy. One 
can easily see that I is symmetric and positive definite operator and consequently, 
one can choose an orthogonal basis in which the operator has the diagonal form 
I = diag(/i ,l2,h). Then I\ , I2 , I3 are called the principal moments of inertia, with 
respect to the principal axes of inertia. If some of Ii,h,h coincide, for example 
if I\ = I2, then any axis in the coordinate plane XOY is principal. The ellipsoid 
(Ifl, fi) = 1 is called inertia ellipsoid of the body at a point O. In the principal 
coordinates its equation is: 



Any symmetry of the body gives the symmetries for the inertia ellipsoid. For 
example the regular hexagon with homogeneous mass distribution is invariant under 
rotations by 7r/3 around the normal line throw the center. Consequently, I\ = I2, 
and any axis in the plane of the hexagon through its center is the principal axis. 
The similar conclusion can be derived for the star (see picture 2). So, we have here 
two geometrically different objects with the same inertia momenta. 



I(u) = \ mid; 



!• =^mj(«x Qi,u x Qi) = I '^m i Qi x (u x Qi),u) = (Iu, u), 




hQl + I 2 n 2 2 +/3O3 = 1. 




Picture 2: The regular hexagon and the star have I\ = I2 



The kinetic energy of the body is given by: 
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Similarly, for the angular momentum M with respect to the point O, we have: 
M = Qj x (rriiVi) = m i Q l x (f2 x Q t ) = 

We consider a motion of a heavy rigid body fixed at a point. Let us denote by 
X the radius vector of the center of masses of the body multiplied with mass m of 
the body and gravitational acceleration g. By T we denote the unit vertical vector. 

The motion in the moving reference frame is described by the Euler-Poisson 
equations [551 H"2"] : 

M = M xCl + f xy 
(3.2) % _ 

r = rxa 

Using that M — I ft, one see that (|3.2[) as a system of six ordinary differential 
equations in M and T with six parameters / = diag(7i, 1%, 13), x — Zq). 
These equations have three first integrals: 

( 3 - 3 ) F 1 = (M,f), 

F 2 = (?,?). 

Since the equations preserve the standard measure, by Jacobi theorem (see for 
example [331 [7]) for integrability in quadratures one needs one more additional 
functionally independent first integral. 

On the other hand, the equations (13.21) are Hamiltonian on the Lie algebra e(3) 
with the standard Lie-Poisson structure: 

(3.4) {M i ,M j } = -e ijk M k , {M h Tj} = -e ijk T k , i,j,k = 1,2,3. 

The structure (|3.4[) has two Casimir functions F\ and F2 from (13.3[) . So, symplectic 
leafs are four-dimensional (they are diffeomorphic to the cotangent bundle of the 
two-dimensional sphere [10]) and for the integrability in Liouville sense one needs, 
besides the Hamiltonian H from p. 31) . one more functionally independent first 
integral. 

From the facts given above, one concludes that a natural problem arises: for 
which values of the parameters I\, I2, h, X , Y 0} Z , the equations (|3.2I) admit the 
fourth functionally independent first integral? 

3.1. Itltegrable cases. Existence of additional independent fourth integral gives 
strong restrictions on moments of inertia and vector x- Such integral exists only in 
three cases: 

• Euler case (1758): Xq = Yq = Zq = 0. The additional integral is F4 = 
(M, M). 

• Lagrange case (1788): I\ = I2, X = (0, 0, Z$). The additional integral is 
F 4 = M 3 . 

• Kovalewski case (1889): h = h = 2I 3 , x = (X ,0,0). The additional 
integral is F 4 = (Of - Q| + jfTx) 2 + (2^i0 2 + jfT 2 ) 2 
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Wc have also cases that admits a fourth integral only with a fixed value of one 
of the integrals. If Casimir function F\ = 0, then we have 

• Goryachev-Chaplygin case (1900): h = h = 4J 3 , X = PQ),0,0). The 
additional integral is F 4 = M 3 (Mf + Mf ) + 2Mir 3 ; 

Beside the completely integrable cases, there are cases that instead of additional 
first integral have an invariant relation. We will focus on Hess-Appel'rot case. Hess 
in [34] and Appel'rot in [5] found that if the inertia momenta and the radius vector 
of the centre of masses satisfy the conditions 



(3.5) Y = 0, 

then the surface 
(3.6) 



X^h{I 2 - h) + Zoy/hih - h) = 0, 



F 4 = MiX + M 3 Z = 



is invariant. It means that if at the initial moment one has that F4 — 0, then this 
will be satisfied during the time evolution of the system. 



3.2. Classical integration of Hess-Appel'rot case. Classical integration of the 

Hess-Appel'rot system is done in the so-called Hess coordinates. The details and 
historical notes can be found in [33] . Hess introduced new coordinates v, p, p: 



Ml 



Mf 



Ml 



(3.7) 



/' : 

9 : 



-(Xorx + Yora + Zors) 

MiXo + M 2 Y + M 3 Z . 



h 
2 



1 /Ml 

2 V7T 



Ml 



Ml 



If one denotes 



Ml 



M, 



Mi 



-Xn 



II 

Mi 
a = — 
h 

Hi = 2p, S 2 
the equations of motion become (see [3"3"]): 

S 2 



1 1 1,1 3 

J 9 Jo 



M 2 
h 



-Yn 



M 3 



y 2 



yl 



1 dv 
2~dl 



1 



S 2 

P 
a 



P 
v 

Pi 



Ml 



H 2 - p 2 )% = (^S 2 - pay-"- + ( Cl 6 2 - ppY-T-, 



p 

P 1 Cl 
P Cl v 

1 dv , r2 s dp 

dt — - ) 2dt + (ciS ~ W) Tf 

where h and ci are fixed values of the first integrals Q3.3p : Fx = c\, H = h. 
The equation of the invariant surface becomes p — 0, and we get: 

S 2 p 
P 1 Cl 
ci v 



1 dv 
2~dl 



dp 
'~dt 



Pi 



1 dv 
2~dl' 
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From the second equation one has /i = cv, where c is a constant, and from the first 
we have 

—1^^^=^^^= = 2dt. 

\J —c 2 6 2 + S 2 v — c 2 v 3 

Hence, v is v = <f>(£ + c), where $(t) is an elliptic function. 

So, in the Hess-Appel'rot case, in the Hess coordinates one can find solutions: 

fj, = c$(t + c), v = $(t + c), p = 0. 

Nevertheless, when Hess-Appel'rot conditions are satisfied, the coordinate trans- 
formation (|3.7p becomes degenerate. Namely, from the system: 

M\ + Ml + Mf = v 

M? M 2 Mn 

Ji H H 

M\Xl - m\zI = 

one needs to find M 2 ,M 2 ,Ml as functions of fi± and v. The determinant of the 
system is: 



A = 



1 



-{Xlh{I 2 -I s )- zlh{h - I 2 )) = 0, 



hhh 

hence, the coordinate transformation (|3.7I) becomes degenerative. Consequently, 
to find a solution of the system one needs to solve one more differential equation. 
Hess proved in 34 that additional differential equations can be reduced to a Riccati 
equation. Let us present Hess result. 

From conditions (|3.6[) . one gets M\ = Z Q u, M 3 = — X u. Let v is defined with 



M 2 = Sv, S 



Zq. From the first integrals 



Ti+T 2 2 +Ti = l, 
XqTi + Z0T3 = — - 



Mf 



M 2 



Ml 



h 
2 



MiFi + M 2 F 2 + Af 3 r 3 = ci = 
by expressing T\, T 2 , T3 and putting them in the first two equations of motion 



I/, = (i-iW 2 M 3 - 
-J3 i-i' 

u 2 = [^--^]M X M 3 



r 2 z 

— TiZq 



one gets 



pZ 



idv + vdu 



~Z„WH 



-vdu — udv 



1 5{u 2 



XoZ 
H 1= S 2 



1 

h 



1 

h 



u 2 du — 5 2 Xq 



(u 2 + v 2 ) 



S 2 - 



1 1 

h~T 2 
(h - V 



uvdv = 0, 



h(h-h) 



X^(u 2 +v 2 ) - h 
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The polynomial Hi is of degree three in p\ — u 2 + v 2 , Hi = P3(p 2 ). Introducing 
u = pi cos ip, v = pi sin tp, the previous equation reduces to 



Lp\ cos ip J dpi = 0, L = 2Xq^o- 1 



2p! ' ri V A/s 

This is the equation derived by Hess. Nekrasov proved that this equation can 
be reduced to a second order linear differential equation with double-periodical 
coefficients. Introducing r = tg%, the last equation becomes a Riccati equation: 

where dilj — s f ldpi , Introducing r = ds / d ^ W e have 

= ^(piW) 

which is the equation obtained by Nekrasov. 

3.3. Lax representation for the classical Hess-Appel'rot system. Algebro- 

geometric integration. The Lax representation for classical Hess-Appel'rot sys- 
tem, with the algebro-geometric integration procedure was presented in [T8 . It is 
proved there that the integration also leads to an elliptic function and an additional 
Riccati differential equation. 

Using isomorphism (|3.ip . equations (13.2[) can be written in the matrix form: 

m = [m, n] + [r, x) 
r = [r,o], 

where the skew-symmetric matrices represent vectors denoted by the same letter. 
We have the following: 

Theorem 3.2. |18) If condition (|3.6I) is satisfied, the equations of Hess-Appel'rot 
case can be written in the form: 



(3.8) , , %2 



L(X) = [L(X),A(X)], 
L(A) = \ 2 C + \M+T, A(\) = x x + n, C = h X - 

The spectral curve, defined by: 

C : p(p, A) := dct(L(A) - pE) = 0, 

is: 

C : -p{p 2 - uj 2 + 2 AA*) =0 
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where 

Xq Zq 

a = — p = — 

V x a + Z o V x a+ Z a 
A = y + Xx, A* = y + Ax, 

(3J) y = io8ri-ar 3 -ir 2 ), i = -^G0Mi - «m 3 - *Af 2 ), 

w = -i [a(dA 2 + Mi A + r x ) + /?(C 3 A 2 + M 3 A + T 3 )] 
= -t [a^A 2 + T x ) + /3(C 3 X 2 + T 3 )] • 

This curve is reducible. It consists of two components: the rational curve C\ given 
by \i — 0, and the elliptic curve Ci : 

(3.10) fj 2 = P 4 (A) = lo 2 - 2AA*. 

The coefficients of the spectral polynomial are integrals of motion. If one rewrites 
the equation of the spectral curve in the form: 

p(fi, A) = -^(/i 2 + AX 4 + BX 3 + DX 2 +EX + F) = 0, 



one gets: 



72 \ 



A = Ii(X 2 + Z 2 j, 
B = 2I 2 {M 1 X + M 3 Z )(= 0), 
D = M 2 + M 2 + M 2 + 2I 2 {X T 1 + z r 3 ), 
E = 2(M 1 T 1 +M 2 T 2 + M 3 r 3 ), 

^ = r 2 + r 2 + r 2 (= 1). 

So, L-A pair (|3.8p gives three integrals and one invariant relation. 

Here, we will review some basic steps in the algebro-geometric integration proc- 
dure from [Tg] , 

Let (/i, f 2 , f 3 ) T denote an eigenvector of the matrix L(X), which corresponds to 
the eigenvalue /i. Fix normalizing condition fx = 1. 
Then one can prove: 

Lemma 3.1. TE\ The divisors of f 2 and f 3 on C 2 are: 

(/ 2 ) = -Pi + P2 - v + 9, 
(/ 3 ) = Pi - P 2 + v - v, 

where P\ and P 2 are points on C 2 over X = 00, and v G C2 is defined with v\ = — 
Vp = -u \x=-i- 

Now we are going to analyze the converse problem. Suppose the evolution in 
time of the point v is known. For reconstructing the matrix L(X), one needs 
x = \x\e larsx , y — \y\e laxsv as functions of time. 

Lemma 3.2. [IS] The point v € T 2 and the initial conditions for M and T deter- 
mine \x\, \y\ and argy — argx, where x and y are given by (|3.9|) . 



16 



BORISLAV GAJIC 



So, in order to determine L(X) as a function of time, one needs to find the 
evolution of the point v and arg a; as a function of time. In |18) the following two 
theorems are proved: 

Theorem 3.3. [18) The integration of the motion of the point v reduces to the 
inversion of the elliptical integral 

dX 1 

— t. 



, ro y/n 2 - 2AA* h 
Denote by <j> x — argx, and u — tg^-. 
Theorem 3.4. |18j The function u(t) satisfies the Riccati equation: 



^ = [f(t) + 9(t)]u 2 + [f(t)-g(t)} 1 



where 



\x\ is a known function of time. 

As we presented in the previous subsection, the classical integration procedure 
also yields one elliptic integral and the Riccati equation. 

3.4. Zhukovski's geometric interpretation. In 47 Zhukovski gave a geometric 

interpretation of the Hess-Appel'rot conditions. Denote Ji — 
Let us consider the so-called gyroscopic inertia ellipsoid: 

M? Ml M\ 
— - + — + — = 1, 
Ji Ji J3 

and the plane containing the middle axis and intersecting the ellipsoid at a circle. 
Denote by I the normal to the plane, which passes through the fixed point O. Then 
the condition (|3.5[) means that the center of masses lies on the line I. 

If we choose a basis of moving frame such that the third axis is /, the second 
one is directed along the middle axis of the ellipsoid, and the first one is chosen 
according to the orientation of the orthogonal frame, then (see [12]), the invariant 
relation (13.61) becomes 



F 4 = M 3 = 0, 

the matrix J obtains the form: 

/ Ji J 13 
J= Ji 

\Jl3 J 3 

and X = (0,0,Z ). 
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One can see here that the Hess-Appel'rot system can be regarded as a perturba- 
tion of the Lagrange top. In new coordinates the Hamiltonian of the Hess-Appel'rot 
system becomes 

Hha = \{Ji{Mf + Ml) + J 3 A/ 3 2 ) + Z T 3 + J 13 M X M 3 = H L + J 13 MiM 3 

This was used as a motivation for a definition of the higher-dimensional Hess- 
Appel'rot systems in [20] . 



4. Higher-dimensional generalization 

Now we will pass to the higher-dimensional rigid body motion. Let us consider 
motion of N points in R™ such that the distance between each two of them is 
constant in time. As an analogy with the three-dimensional case, we have two 
reference frames: the fixed and the moving ones. In the moving reference frame, 
the velocity of the z-th point is: 

Vi(t) = B- l qi{t) = B~ x BQi = n{t)Q, 

where again Qi represents the radius vector of the i-th point, and £1 is skew- 
symmetric matrix (fl £ so(n)) representing the angular velocity of the body in 
the moving reference frame. The angular momentum is a skew-symmetric matrix 
defined by 

i i 

= 53 m l {nQ l Q t l + QiQln) = ni + m, 

i 

where / = ^\ rriiQiQl is a constant symmetric matrix called the mass tensor of 
the body (see [50]). 

If one chooses the basis in which / = diag(/i, ...,/„), the coordinates of angular 
momentum are My = (I; + 
The kinetic energy is 

i i 

Since it is a homogeneous quadratic form of angular velocity f2, one has (|^, fi) 
27 where (A, B) = -\Trace{AB) is an invariant scalar product on so(n). One 
gets 

grp 

m 

or |^ = M and finally 
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Remark 4.1. The expression My = (7j + gives a left-invariant metric on 

SO(n). In that sense, the solutions of the classical Euler equations can be inter- 
preted as geodesic lines of the left-invariant metric on 50(3). Arnol'd generalized 
the Euler equations (see [6]). He derived equations of geodesies of an arbitrary 
left-invariant metric on Lie group G. In Manakov found L-A pair for the wider 
class of metrics Mij — flij, and showed that this class belongs to the class 

considered by Dubrovin in [25] , and hence, the solutions can be expressed in theta 
functions. 

The Lie group E(3) can be regarded as a semidirect product of the Lie groups 
SO{3) and K 3 . The product in the group given by 

(A lt n) ■ (A 2 ,r 2 ) = {A 1 A 2 ,r 1 +A 1 r 2 ) 

corresponds to the composition of two isometric transformations of the Euclidian 
space. The Lie algebra e(3) is a semidirect product of R 3 and so(3). Using iso- 
morphism between the Lie algebras so(3) and R 3 , given by (I3.1[) . one concludes 
that e(3) is also isomorphic to the semidirect product s = so(3) x a d so(3). The 
commutator in s is given by: 

[(ai,6i), (03,62)] = ([01,02], [01,62] + [61, o 2 ]). 

One concludes, that there are two natural higher-dimensional generalizations of 
equations (|3.2p . The first one is on the Lie algebra e(n) that is a semidirect product 
of so(n) and R n . The n-dimensional Lagrange case on e(n) is defined in |i9 , where 
its integrability is proved. The higher-dimensional Kowalevski case together with 
Lax representation is constructed in [TT] (see also [5]). For a list of integrable cases 
see for example |46j . 

The second one, given by Ratiu in [44] is on semidirect product s = so(n) x ac ; 
so(n). Equations of motion in moving frame are ( |44| ) : 

(4.1) M = [M,n] + [r, x ], r = [r,fi]. 

Here M £E so(n) is the angular momentum, Q 6 so(n) is the angular velocity, 
X G so(n) is a given constant matrix (describing a generalized center of the mass), 
F G so(n) . Angular momentum M and f2 are connected by M = Ifl + OJ. If the 
matrix I is diagonal, I — diag(/i, . . . , J n ), then My = (Jj + ■ The Lie algebra 

s is the Lie algebra of Lie group S — SO(n) XAd so(n) that is semidirect product of 
SO(n) and so(n) (here so{n) is considered as the Abelian Lie group). The group 
product in S is (Ai,h) ■ {A 2 ,b 2 ) = (A 1 A 2 ,b 1 + Ad Al b 2 ). 

Ratiu proved that equations (|4.ip are Hamiltonian in the Lie-Poisson structure 
on coadjoint orbits of group S given by: 

{/) ff}(M) v )=~ Ki d if(Vi v),dig{n, 

(4.2) -v(WfM&gfav)]) 

- v{[d2f{n,v),d\g{ii,v)]), 

where f,g are restrictions of functions / and g on orbits of coadjoint action and 
dif are partial derivatives od df. On so{n) a bilinear symmetric nondegenerate 
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biinvariant (i.e. &([£,??],£) = fc(£, [17, C] ) ) two form exist, which can be extended to 
s as well: 

k a (($i,vi), (£2,172)) = H£,i,m) + H&iVi)- 

Hence, one can identify s* and s. Then, the Poisson structure (|4.2p can be written 
in the form 

{/,§}(£. V)=~ [(ffroda/)^, (Wi<?)(£, *?)]) 

(4-3) - fc(£, [(5rodi/)fo r,), (5rad 25 )(e,7?)]) 

- fc(?7, [(ffracfe/Xf , f), (grad 2 g)(£,, 11))), 

where gradi are fc-gradients in respect to the i-th coordinate. 

In [44], the Lagrange case was defined by I\ = I 2 = a, I3 = ••• = I n = 
b, X12 = -X21 f 0, Xij = 0, i {(1,2), (2,1)}. The completely symmetric 

case was defined there by I\ = ■ ■ ■ = I n = a, where \ € so(n) is an arbitrary 
constant matrix. It was shown in [33] that equations (|4.ip in these cases could be 
represented by the following L-A pair: 

^(A 2 C + AM + r) = [A 2 C + AM + r, A% + f2], 

where in the Lagrange case C = (a + b)x, and in the symmetric case C = 2a%. 

4.1. Four-dimensional rigid body motion. To any 3x3 skew-symmetric matrix 
one assigns one vector in three-dimensional space using isomorphism between R 3 
and so(3). Using the the isomorphism between so(4) and so(3) x so(3), one can 
assign two three-dimensional vectors A\ and A2 to (4 x 4)-skew-symmetric matrix 
A. 

Vectors Ai and A2 are defined by: 



Ai = 



A+ + A- 



A 2 = 



A+-A. 



where A+, A- € 



(4.4) 



2 ' * 2 
I 3 correspond to Aij £ so(4) according to: 



( 


-A\ 













-A 2 









-A 3 




A 2 . 


A 3 . 


/ 



Here A± are the j-th coordinates of the vector A± . 

By direct calculations, we check that vectors 2A\ x B\ and 2A 2 x f? 2 corre- 
spond to commutator [A, B], if vectors Ai,A 2 and Bi, B 2 correspond to A and B 
respectively. 

Consequently, equations of motion (|4.1[) on so(4) x so(4) can be written as: 

Mi = 2(Mx xfii+Lx Xl ) f ! = 2(Px x n x ) 

M 2 = 2(M 2 x ft 2 + L 2 x X2 ) T 2 = 2(P 2 x n 2 ) 

Recall that M = Iil + ill. The matrix elements of the mass tensor of the body I 
are Iki = J2 m iQ(i)kQ(i)h k,l — 1,...,4. Choose the coordinates (X±, X 2 , X 3 , X4) 
of the moving reference frame in which / has diagonal form / = diag(ii, I 2 , 13, 1 a). 
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Then, for example h = Yl m i X {i)v 12 = J2 m i X (i) 2 > H m i x {i)i x (i)2 = etc., 
where is the A:-th coordinate of i-th point. In the three-dimensional case the 

moments of inertia were defined with respect to the line through the fixed point 
O. We derive the angular velocity uj as a skew-symmetric matrix the elements of 
which correspond to the rotations in two-dimensional coordinate planes. Hence, 
here it is natural to define the moments of inertia of the body with respect to the 
two-dimensional planes through the fixed point. For example the moment of inertia 
with respect to the plane X\OX 2 is I\ + 1 2 , and M 12 = {h + -^2)^12, etc. 

Here we observe a complete analogy with the three-dimensional case. For ex- 
ample, the moment of inertia with respect to OZ axis I 3S — J2i m i( x f + Y 2 ) 
consists of two addend Y^i m % x 1 an d J2i m iY 2 that are diagonal elements of the 
mass tensor of the body. 



For vectors M + and M_ one has 



Mi 



M_ 



Finally, one can calculate 
Afi 
M 2 



(4.6) 



Z_)fia 



(J+-/_)0 2 ) 
(J++/_)0 2 ) 



At a glance it looks that (|4.5|) are equations of motion of two independent three- 
dimensional rigid bodies. However, the formulas (|4.6p show that they are not 
independent and that each of Mi, M 2 depends on both Vt\ and il 2 . 

4.2. Lagrange bitop. Definition and Lax representation. Generalizing the Lax 

representation of the Hess-Appel'rot system, the new complete integrable four- 
dimensional rigid body system is established in [18| . A detailed classical and 
algebro-geometric integration was presented in [19] . 

The Lagrange bitop is four-dimensional rigid body system defined by (see [181 

my- 



(4.7) 



h = I 2 = a 
h = h = b 



and 








X12 





°\ 




-X12 






















X34 










-X34 


J 



with the conditions a ^ b, X12, X34 0, |xi 2 | 7^ |X34|- 
We have the following proposition: 

Proposition 4.1. p75l IT9"] The equations of motion (14. ip under conditions (|4.7|) 
have an L — A pair representation L(X) — [L(X), A(X)] , where 



(4.8) 

and C = (a 



l(x) = x 2 c + am + r, A(x) = x x + n, 



b) X - 
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Let us briefly analyze spectral properties of the matrices L(X). The spectral 
polynomial p{X 1 fj,) = det (L(X) — fj, ■ 1) has the form 

j>(A,aO =n i + P(X) f i 2 + lQ(X)} 2 , 



where 
(4.9) 



P(X) = AX i + BX 3 + DX 2 + EX + F, 
Q(A) = GA 4 + HX 3 + IX 2 + JX + K. 



Their coefficients 

A = C\ 2 + G 3 2 4 = (G+, G+> + (CL, C-), 

B = 2G 34 Af 34 + 2C 12 M 12 = 2 ((C+,M+) + (C_,M_)) , 

£> = M 2 3 + M 2 4 + M| 3 + Afg, + M| 4 + 2Gi 2 r 12 + 2G 3 4r 3 4 

= (M+, M+) + (M_, M_) + 2 ((G+, r+) + (c_,r_)) , 
£ = 2r 12 Afi 2 + 2ri 3 Afi 3 + 2r 14 Afi4 + 2r 23 M23 + 2r 24 Af 24 + 2r 34 M 3 4 

= 2((r + ,Af + > + (r_ ; Af_)), 
^ = r? 2 + r 2 3 + r 2 4 + r 2 3 + r 2 4 + r 2 4 = (r+, r+) + (r_, r_), 

G = Gi2G 34 = (G+,G_), 

H = C 3i M u + C 12 M M = (G+, M_) + (G_,M+), 
7 = G 34 r 12 + r 34 Gi2 + Mi 2 M 34 + M 23 Mi 4 - M 13 M 2i 

= (c+,r_) + (G_,r+> + (m+,m_), 
j = Af 34 r 12 + M 12 r 34 + Mi 4 r 23 + M 23 r 14 - r 13 M 24 - r 24 Af 13 

= <M + ,r_) + (M_,r + ), 
if = r 34 ri2 + r 23 Ti 4 — ri 3 r 24 = (r+,r_). 

are integrals of motion of the system (|4.1I) , (|4.7[) . Here M + , M _ S K 3 are defined 
with (|4.4I) (similar for other vectors). System (|4.f p . (|4.7[) is Hamiltonian with the 
Hamiltonian function 

H = -(Miafiia + Mi 4 14 + Af 23 ft 23 + Af 12 ft 12 + A/ 34 tt 34 ) + xi 2 ri 2 + Y 34 r 34 . 

The algebra so(4) x so(4) is 12-dimensional. The general orbits of the coadjoint 
action are 8-dimensional. According to [44], the Casimir functions are coefficients 
of A , A, A 4 in the polynomials [det L(A)] 1 / 2 and — \Tr(L(X)) 2 . One calculates: 

detL(X) V2 = GA 4 + HX 3 + IX 2 + JX + K, - ^Tr (i(A)) 2 = AX 4 + EX + F. 

So, Casimir functions are J, K, E, F. Nontrivial integrals of motion are B, D, H, I, 
and, one can check that they are in involution. When |xi 2 | = | X34 1 , then 2H = B 
or 2H = —B and there are only 3 independent integrals in involution. Thus, 

Proposition 4.2. [H] For \xiz\ 7^ \X3i\i system (|4.f [) . (|4.7I) is completely integrable 
in the Liouville sense. 
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System (|4.1I) , (|4.7|l doesn't fall in any of the families defined by Ratiu in [33] 
and together with them it makes complete list of systems with the L operator of 
the form 

L(\) = X 2 C + XM + T. 

More precisely, if X12 ^ 0, then the Euler-Poisson equations (|4.ip could be written 
in the form (I4.8[) (with arbitrary C) if and only if equations (|4.1I) describe the 
generalized symmetric case, the generalized Lagrange case or the Lagrange bitop, 
including the case X12 = ±X34 [IB] , 



4.2.1. Classical integration. For classical integration we will use equations (|4.5[) . 
On can calculate that 

Xi = (0, 0, --(xi2 + X34)), X2 = (0, 0, -~(xu - X34)) 

and also 

Mi = ((a + 6)ft (1)1 , (a + 6)0 (1)2 , (a + 6)fi (1)3 + (a - o)Q (2)3 ) 
M 2 = ((a + 6)0 (2)1 , (a + &)fi (2 ) 2 , (a - &)%)3 + (a + ^(2)3)- 

If we denote S7i = (pi, (ft, Jl), £7 2 = (p 2 , g 2 , r 2 ), then the first group of the 
equations f|4.5[) becomes 

pa - m<7ir 2 = -niF (1 ) 2 , p 2 - mq 2 ri = -n 2 F( 2 ) 2 

<ji + mpir 2 = niT (y!, q 2 + mp 2 ri = n 2 F (2)1 

(a + b)fi + (a - b)r 2 =0, (a - b)f x + (a + o)r 2 = 

where 

2(o - 0) 2x(i) 3 2x (2 )3 

m = T7T' ni = TT' 712 = ~T' 

a + o a + o a + o 

The integrals of motion are for i = 1,2: 
(a + b)aiX(i)3 = fn 

(a + 6)[(o + b){pf + qf) + (a + b)aj + 2x (l )3r (l)3 ] = fa 
(a + b)piT^i + (a + 6)g,T(i) 2 + (a + 6)a 4 r W3 = / i3 

r (i)i + r («)2 + r (j)3 = 1) 



where 



(a + 6)ri + (a — 6)r 2 (a + 6)r 2 + (a — 6)n 

ai = rr a 2 = -t 

a+o a+o 

(a + 
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Introducing Pi,<Ji, defined with pi — pi cos <7j, qi = piSmat, after calculations, one 
gets 

/l3 



2-i 2 / ./ 13 \ 

p x ox + mr 2 p 1 = nii——^- ~ air (1)3 ) 



(4.10) l( P 2 y} 2 =4n 2 P 2 {l-\(a l+ p 2 ) 2 } -4^(4^7 ~a iai -%lf, i = 1, 2 



2 • i 2 / /23 \ 

p 2 a 2 + mri/? 2 = n 2 — — - a 2 r ( 2) 3 )- 

Let us denote Ui = p\ , u 2 = p 2 .- From (|4. 10[) we have 

u 2 l =P l {u i ), i = 1,2, 
P;(u) = -4u 3 - 4w 2 B. i + 4uCi + Di, i= 1, 2; 

Bi = 2ai + af, d = nj - af - 4 y - 2afa 

(a + o) z 

= ~ 4 (~7 — TTTT + a * a i) ' 1 = 1,2. 



(a + &) 2 

From the previous relations, we have 



dm r du2_ = L 



So, the integration of the Lagrange bitop leads to the functions associated with the 
elliptic curves E 1 ,E 2 where E t = Ei{au <ij, X(i)3, fi2, fa) are given with: 

(4.11) E z :y 2 = P t {u). 

Equations f)4.5[) are very similar to those for the classical Lagrange system. How- 
ever, the system doesn't split on two independent Lagrangian systems 

4.2.2. Properties of spectral curve. The spectral curve is given by: 

C: p 4 + P(\)p 2 + [Q(\)} 2 = 



where P and Q are given by (J4T9J) . 

There is an involution a : (A, fi) — > (A, —fi) on the spectral curve which cor- 
responds to the skew symmetry of the matrix L(X). Denote the factor-curve by 
d=C/a. 

Lemma 4.1. [19] 

• The curve C\ is a smooth hyperelliptic curve of the genus g(Ci) = 3. The 
spectral curve C is a double covering of C± . The arithmetic genus of C is 

9a(C)=9. 

• The spectral curve C has four ordinary double points Si,i = 1, . . . , 4. The 
genus of its normalization C is five. 

• The singular points Si of the curve C are fixed points of the involution a. 
The involution er exchanges the two branches of C at Si. 
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In general, whenever matrix L(A) is skew-symmetric, the spectral curve is re- 
ducible in an odd-dimensional case and singular in an even-dimensional case. 

The detailed algebro-geometric integration procedure of the system is given in 
|19j . Analysis of the spectral curve and the Baker- Akhiezer function shows that 
the dynamics of the system is related to a certain Prym variety II that corresponds 
to the double covering defined by the involution a and to evolution of divisors of 
some meromorphic differentials fij. It appears that 



<4 



are holomorphic during the whole evolution. Compatibility of this requirement 
with the dynamics puts a strong constraint on the spectral curve: its theta divisor 
should contain some torus. In the case presented here such a constraint appears 
to be satisfied according to Mumford's relation. These conditions create a new 
situation from the point of view of the existing integration techniques. For details 

see nm. 



4.3. Four-dimensional Hess-Appel'rot systems. The starting point for construc- 
tion of generalization of the Hess-Appel'rot system was Zhukovski's geometric inter- 
pretation given in subsection !3.4l Having it in mind, in |20) the higher-dimensional 
Hess-Appel'rot systems are defined. First we will consider the four-dimensional 
case on so(4) x so(4). We will consider metric given with ft — JM + MJ. 

Definition 4.1. [50] The four-dimensional Hess-Appel'rot system is described by 
equations (|4.1[) and satisfies the conditions: 

(1) 



(4.12) 



O = MJ + JM, J = 



Jl 





Jl3 


M 





Jl 





J24 


Jl3 





■h 





V° 


J24 








(2) 



( 





X12 





°\ 




X12 






















X34 


\ 








-X34 


0/ 



The invariant surfaces arc determined in the following lemma. 

Lemma 4.2. [20 For the four- dimensional Hess-Appel'rot system, the following 
relations take place: 

M12 = Ji 3 (M 13 M 12 + M24M34) + J24(M 13 M 34 + Mi 2 M 24 ), 

M 34 = Ji 3 (-Mi 3 M34 - M 12 M 24 ) + J2i(-M 13 M 12 - M 24 M 3i ). 

In particular, if Myi = M34 = hold at the initial moment, then the same relations 
are satisfied during the evolution in time. 
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Thus, in the four-dimensional Hess-Appel'rot case, there are two invariant rela- 
tions 

(4.13) Mia = 0, M 34 = 0. 

Let us now present another definition of the four-dimensional Hess-Appel'rot 
conditions, starting from a basis where the matrix J is diagonal in. 
Let J = diag( Ji, J 2 , J3, J4). 

Definition 4.2. [20] The four-dimensional Hess-Appel'rot system is described by 
the equations (|4.1I) and satisfies the conditions: 

(1) 

il = MJ + JM, J = diag(Ji, J 2 , J 3 , J 4 ), 

(2) 



(3) 



where 






X12 







Xl2 





X23 








-X23 





X34 


Xl4 





-X34 


/ 




• J 4 = 


■^2 — Jl, 




J 3 - 


Jl 


Ji — J 2 






-*? 


Vi + *2 





h := 



i 2 



2(X14X34 - X12X23) 



Xl4 



X34 



Xl2 



~2 ' 

X23 



2(Xi4Xi 



X23X34J 



"X?4 



X§4 



Xl2 



x| 3 ' 



Proposition 4.3. [20] There exists a hi- correspondence between sets of data from 
Definition ^. 1\ and Definition ^. 2\ 

Remark 4.2. 1) In the case J24 7^ 0,X34 = 0, there is an additional relation 
X12X34 + X14X23 = 0. It follows from the system 

X12 siny + X23 cos(^ = 0, 

Xi 4 sin<y9 - X34 cos(^ = 0, 

2) In the case J24 = 0,^34 = 0, additional relations are X34 = Xi4 = 0, an d the 
second relation from Definition 14.21 can be replaced by the relation 



X12V J2- Ji+ X23 1 

Theorem 4.1. |20j The four- dimensional Hess-Appel'rot system has the following 
Lax representation 

L(X) = [L(X),A(X)], 

1 



Ja - J 2 



0. 



L(A) 



A 2 C + AM + T, A(A) = Ax + fi, C 



Ji + J; 



2<i 



BORISLAV GAJIC 



where 



One can calculate the spectral polynomial for the four-dimensional Hess-Appel'rot 
system: 

p(A, p) = det(L(A) - fi • 1) = u 4 + P(X)p 2 + Q(A) 2 , 

P(A) = aA 4 + 6A 3 + cA 2 + dX + e 
Q(A) = /A 4 + .gA 3 + h\ 2 + iX + j 

a = Cl2 + ^34; 

6 = 2C 12 M 12 + 2C 34 M34(= 0), 

c = Af 2 3 + Af 2 4 + M 2 3 + M 2 4 + M 2 2 + M 2 4 + 2C 12 T 12 + 2C 3 iT 3i , 

d = 2r 12 M 12 + 2r 13 Mi 3 + 2r 14 Mi4 + 2r 23 M 23 + 2r 24 M 24 + 2r 34 M 34 
e = r 42 + rj 3 + r 2 4 + r 2 3 + r 2 4 + r 2 4 , 

/ = Cl 2 C 3 4 

9 = C 12 M 3A + C 34 Mi 2 (= 0), 

ft = r 34 c"i 2 + ri 2 c 34 + M12M34 + M 23 Mi4 - M13M24, 

i = M 34 ri 2 + Mi 2 r 34 + M M r 23 + A/ 23 r M - ri 3 M 24 - r 24 Mi 3 , 

j = r 3 4ri 2 + r 23 ri4 — Ti 3 r 2 4- 

In the standard Poisson structure on semidirect product so(4) x so(4) the functions 
d,e,i,j are Casimir functions, c, h are first integrals, and b = 0, g = are the 
invariant relations. As we already mentioned general orbits of co-adjoint action 
are eight-dimensional, thus for complete integrability one needs four independent 
integrals in involution. 



4.4. The n-dimensional Hess-Appel'rot systems. In 20] we introduced also 

Hess-Appel'rot systems of arbitrary dimension. 

Definition 4.3. The n-dimensional Hess-Appel'rot system is described by the 
equations (|4.ip . and satisfies the conditions: 

(1) 



n = jm + mj, j = 



fJl 





Jl3 





. 







Jl 





J24 


. 


. 


Jl3 





■h 





. 


. 





J 2 4 





h 


. 


. 














. 


. 


u 











. 


. J 3 ) 
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(2) 



/ 





Xl2 


. 


• °\ 




-Xl2 





. 


. 










. 


. 










. 


. 


\ 








. 





Direct calculations give the following lemma: 

Lemma 4.3. |20) For the n- dimensional Hess-Appel'rot system, the following re- 
lations are satisfied: 

(1) 

n 

Mia = Ji3(M 12 M 13 + M24M34 + M 2 pM 3p )+ 

p—5 



J 24 (M 12 M 24 + M 13 M 34 - Y M lp M, 



-ip) 



p—5 



M: 



34 



- Ji 3 (M 13 M 3 4 + M 2i M 12 + M lp M pi ) 

p—5 



J 24 (M 13 M 12 + M24M34 + J2 M 2 P M 3p ), 

p— 5 

M 3p = -Ji 3 (M 13 M 3p + M 2p M 12 ) - J 24 (M 34 M 2p + M 23 M 4p ) + 



M u n 4p - n 34 M 4p + J2{M 3k n kp - n 3k M 4p ), P >4, 

Mi P = Ji3(-M u M 3p + M lp M 34 ) + J 24 {M 12 M lp - M 24 M 4p )- 

n 

M M n 3p + ci m m 3p + ^{M 4k n kp - n 4k M 4p ), P > 4, 



k=5 



(2) 



M M = 0, k,l>4. 

(3) The n-dimensional Hess-Appel'rot case has the following system of invari- 
ant relations 

M 12 =0, M ip = 0, l,p>3. 

By diagonalizing the matrix J, we come to another definition 

Definition 4.4. [5D] The n-dimensional Hess-Appel'rot system is described by the 
equations (|4.ip . and satisfies the conditions 

(1) 

fi = JM + MJ, J = diag(Ji, J 2 , J 3 , J 4 , J 4 ), 



2,x 
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(2) 



(3) 






Xl2 





Xl4 • 


■ °\ 


-Xl2 





X23 





. 





"X23 





X34 • 


. 


~Xl4 





~X34 


. 


. 











. 


• V 


J 3 


- J 4 


= J 2 - 


Jl, 






- Jl 










+t\ 


x/i- 


vt\ 





where 



t 2 



X12X34 + X14X23 = 
2(Xl4X34 ~ Xl2X2 3 ) 

Xu ~ X34 + X?2 - X23 ' 

2(Xl4Xl2 - X23X34) 
-Xu - X§4 + X?2 + xh ' 



As in the dimension four, there is an equivalence of the definitions. 

Proposition 4.4. [20] There exists a bi- correspondence between sets of data from 
Definition \4.3\ and Definition ^ -4\ 

The following theorem gives a Lax pair for the n-dimensional Hess-Appel'rot 
system. 

Theorem 4.2. [20] The n- dimensional Hess-Appel'rot system has the following 
Lax pair 

L(X) = [L(\),A(X)}, 

1 



L(A) = YC + AM + T, A(X) = \ x + tt, C 



-X- 



J1 + J3 

4.5. Classical integration of the four-dimensional Hess-Appel'rot system. 

Detailed classical and algebro-geometric integration procedures for the four-dimen- 
sional Hess-Appel'rot case are presented in [20 . Here again equations (|4.5p are 
useful for classical integration. We have: 

Xi = (0, 0, -^(Xi2 + X34)), X2 = (0, 0, --(X12 - X34))- 
Integrals of the motion are 

(M i ,M i )+2—^—( X i,T i ) =hi, 

(4.14) (r l ,r,) = i, i 

(Mi,Ti) = Ci, 
( X i,Mi}=0. 



1,2, 
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Here the metric that gives connections between M and Q is different from that for 
Lagrange bitop. We have 

fti = ((Ji + J 3 )M (1)1 - (J 13 - J 24 )Af (2)3 , (J x + J 3 )M (1)2 , 
(Ji + J 3 )Af(i) 3 + (Ji - J 3 )M (2)3 - (J 13 + J 24 )M (2)1 ), 

^2 = ((Ji + J 3 )Af (2)1 - (Ji 3 + J 24 )M (1)3 , (J x + J 3 )M (2)2 , 
(Ji + J 3 )M (2)3 + (Ji - J 3 )M (1)3 - (J13 - J 24 )M (1)1 ), 

where again Af(j)j is the j-th component of the vector Mj. Using these expressions, 

equations (I4.5[) can be rewritten in the following form: 

(4-15) 

M(i)i =2[(Ji - ^ 3 )Af( 1 ) 2 M( 2)3 - (Ji 3 + J 24 )M( 1 ) 2 Af (2)1 + r(!) 2 X(i)3], 
M (1)2 =2[-(Ji - J 3 )M (2)3 M (1)1 - (Ji 3 - J 24 )A/ (1)3 A/ (2)3 + 

(Ji3 + ^2 4 )Af(i)iM( 2)1 - T( 1 ) 1 X(i) 3 ], 
Af(i) 3 =2(Ji3 - J 24 )Af(i) 2 M( 2 ) 3 , 

f ( i)i =2[r ( i )2 ((Ji + J 3 )Af( 1)a + (Ji - Ja)M m - (Ji3 + J24)M m )~ 

T(i) 3 (^i + >^ 3 )Af(i) 2 ], 
f (i )2 =2[r ( i )3 ((Ji + J 3 )M (1)1 - (J 13 - J 24 )Af (2)3 )- 

r ( i)i((Ji + J 3 )Af (1)3 + (Ji - J 3 )Af (2)3 - (J13 + J 24 )A/ (2)1 )], 

f (i )3 =2[r ( i)i(Ji + j 3 )a/ (1)2 - r ( i )2 ((Ji + j 3 )m (1)1 - (Ji 3 - J 24 )M( 2 ) 3 )], 

and 
(4.16) 

A^(2)i =2[(Ji - J 3 )Af( 2 ) 2 Af(i) 3 - (Ji 3 - J 24 )A/ (2)2 Af (1)1 +r (2)2 X( 2 ) 3 ], 
M(2)2 =2[-(Ji - J 3 )A/ (1)3 A/ (2)1 - (Ji 3 + J 24 )M (2)3 M (1)3 + 
(Ji3 - J 24 )A/( 2 ) 1 Af (1 ) 1 - r( 2 )iX( 2 ) 3 ]) 

M(2)3 =2(Jl3 + J 24 )Af( 2 ) 2 A/(!)3, 

f (2)1 =2[r (2)2 ((Ji + J 3 )A/ (2)3 + (Ji - J 3 )M ( i )3 - (J13 - J 24 )A/ (1)1 )- 

r (2)3 (Ji + J 3 )Af (2)2 ], 
T(2) 2 =2[r (2)3 ((Ji + J 3 )Af (2) i - (J13 + J 24 )Af (1)3 )- 

r (2) i((Ji + J 3 )Af (2)3 + (Ji - J 3 )M C i )3 - (J13 - J 24 )A/ (1)1 )], 
T( 2 ) 3 =2[r (2) i(Ji + J 3 )A/ (2)2 - r (2)2 ((Ji + J 3 )M (2) i - (Ji 3 + J 24 )Af (1)3 )]. 

One can see here that Af( 1 ) 3 = Af( 2 ) 3 = 0, giving two invariant relations introduced 
before. 

Let us introduce coordinates K; and L as follows: 



Af(i)i = Ki sink, M {i)2 = Kicosk, i=l,2. 
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From equations (I4.15[) , (|4.16l) . using integrals (|4.14[) . we have 

2 



r(i )3 = 4(Ji + J 3 ) 2 



(1 - rfo 3 )(fc x 



J1 + J3 



X(l)3 r (l) 3 ) 



Thus r(!) 3 can be solved by an elliptic quadrature. Also from the energy integral 
we have that 

2 



Kf = hi 



Jl+Jc 



■X(i)3r(i)3- 



Since tan l\ 



and 



M (1)2 



we have: 



Zi = -2(Ji 3 + J2a)K 2 sinZ 2 



2X(i)3Ci 



^i 2 



x 2 r 2 1)2 - 2 Cl M (1)2 r (1)2 + c? - M ( 2 in (i - r 2 1)3 ) = 0. 

Similarly, one gets: 

2 r ^ 

"X(2)3 i (2)3j _ 



1(2)3 = 4 W 



i ~t~ ^3) 

2 



(1 



(2)3 



Jl + J3 



= ft(T(a)3), 



Ji + J; 



"X(2)3- 1 (2)3, 



Z 2 = -2(Ji 3 - J24)^i sinZi 



2X(2)3C2 



*2 



^2 r (2)2 - 2c 2 M (2 )2r ( 2)2 + C 2 - M (2)1 (l - I^g) = 0. 

From the previous considerations one concludes that integration of the four- 
dimensional Hess-Appel'rot system leads to a system of two differential equations 
(for l\ and I2) of the first order and two elliptic integrals, associated with elliptic 
curves E± and E% defined by 

Ei'. y 2 =Pi(x) = 8A t x 3 - AB.x 2 - 8A t x - AC t , i = l,2 

where 

Ai = (Ji + J 3 )X(i)3, B % = (Ji + J 3 ) 2 h t , Q = (Ji + J 3 ) 2 (c 2 - hi). 

This is a typical situation for the Hess-Appel'rot systems that additional integra- 
tions are required. 

In [3D] the algebro-geometric integration procedure is presented. It is closely 
related to the integration of the Lagrange bitop. 

In [351 ISS] the Hess-Appel'rot systems are considered within the framework of 
partial reduction. 

5. Another Lax representation for the classical Hess-Appel'rot case. 

Generalizations 

In 1846 Jacobi gave an algebraic description of the Jacobian of a hyperelliptic 
curve. Beauville noticed that using given description, any hyperelliptic curve can 
be seen as a spectral curve of some matrix L(X). Starting from the spectral curve 
for the Lagrange top, given by Ratiu and van Moerbeke, Gavrilov and Zhivkov in 
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[3"2] have constructed a new L-A pair for the Lagrange top. Using a modification of 
that construction, starting from elliptic curve (|3.10p we have constructed another 
Lax representation for the Hess-Appel'rot system. This was a starting point for 
construction of a class of systems on e(3), with the same elliptic curve as a spectral 
curve (for details see [H]). Let u>, A and A* are given by 



Proposition 5.1. On hypersurface (|3.6[) the equations of the Hess-Appel'rot sys- 
tem are equivalent to the Lax representation 



(5.1) 

where 

(5.2) 

and a = a^+pn. 



L(A) 



L(X), 



A 2 L(A) - a 2 L{a) 
A — a 



V2i 



A" (A) 



A(A) 



"A 2- 



Starting from Lax representation (|5.1[) in [3T] the family of new systems is de- 
scribed. Let us consider the general case of equations (|5.ip . with a as an arbitrary 
polynomial in M's and L's. 

The corresponding spectral curve is the elliptic curve 



(5.3) 



w 2 (A) A(A)A*(A) 



A 4 A 4 

Here M ad L are in e(3) as before. In terms of x,y,x,y,x\,y\ from (|3.9p and 
yi = aT\ + PT3, x\ = a.M\ + /3M3, the standard Poisson structure p.4[) on e(3) 
has the form: 

{x,y} = 0, {x,y} = 0, {x, x x } = ix, 
{y, xi} = -iy, {yi,xi} = 0, {x, y%} = 



{x,xi} = -ix, {y,x 1 } = iy, 

-iy, {x,yi}=iy, {yi,v} = 0, 



{yi,y}=0, {x,x} = -ixi, {y,y} = 0, {x,y} = -iy x , {x,y}=iyi. 
Observe that matrices L given by (|5.2D satisfy 



where 



L(A),l(/x) 



L(X) = L(A) 



r(A-/i),i(A) + l(/x) 



1 
1 



ht*) = 



<8L(/i), 



with the permutation matrix as an r-matrix 



KA) = ± 



10 
10 
10 
1 



32 



BORISLAV GAJIC 



Equations (|5.1I) can be rewritten in the form 
Mi = Z T 2 + aZ M 2 , 
M 2 = X T 3 - Z Q T X + a(X M 3 - Z M 1 ), 
M 3 = -XqT 2 — aX M 2 , 

, s • r 2 M 3 - r 3 M 2 

(5.4) ri = 2 3 3 2 + aZo r 2 , 

■ r 3 Mi-riM 3 , 

r 2 = ^ h a(A r 3 - z rij, 

TiM 2 - T 2 M 1 

r 3 = 1 % 2 1 - ax r 2 . 

We have the following Proposition. 
Proposition 5.2. System (|5.4D can oe rewritten as: 
M i = {M i ,H 1 } + a{M i ,H 2 }, 
r i = {T i ,H 1 } + a{T it H a }, 1 = 1,2,3, 



where 



Hi = M i + ^1 + M i + (XoFi + Zo p 3)j ff2 = XqMi + ZoMa 



As we have already mentioned the Poisson bracket (|3.4I) has two Casimir func- 
tions: 

F l = MiTi + M 2 T 2 + A/ 3 r 3 , F 2 = V 2 + T 2 . + T 2 . 

Thus, a symplectic leaf, defined by conditions F\ = ci, F 2 = c 2 is a four-dimensional 
manifold. For integrability in the Liouville sense on e(3), another first integral be- 
side the Hamiltonian is necessary. On the other hand, if a system is not Hamilton- 
ian, generally speaking, five first integrals of motion for integrability in quadratures 
are required. But, if a nonhamiltonian system has an invariant measure, then, ac- 
cording to the Jacobi theorem, for integrability in quadratures one needs four first 
integrals of motion. 

For a general polynomial o, the system (|5.4p is neither Hamiltonian in the Pois- 
son structure (|3.4p , nor preserves the standard measure. A simple criterion for 
preserving the standard measure is given by: 

Proposition 5.3. 21J System (|5.4p preserves the standard measure if and only if 
the polynomial a satisfies the condition: 

{a,X M 1 + Z M 3 } = 0. 
As a consequence we have: 

Proposition 5.4. 21j In the following five cases, the standard measure is preserved 

(i) : if the polynomial a is a Casimir function: a = MiTi + M 2 T 2 + M 3 r 3 ; 

(ii) : if the polynomial a is a Casimir function: a = T 2 + T 2 , + T|; 

(iii) : ifa = X Mi+Z M 3 ; 
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(iv) : ifa = X T 1 +Z r 3 ; 

(v) :ifa= M\ + Mf + M| . 

Theorem 5.1. |21j If Xq ^ 0, or 7^ 0, m the first three cases given above, the 
systems are Hamiltonian, while in the fourth and the fifth cases, the systems are 
not Hamiltonian in the Poisson structure (13.41) . 



If a is a Casimir function, for an arbitrary function / we have 

{/, Hx} + a{f, H 2 } = {/, H X + aH 2 }. 

Hence, in the first two cases the systems are Hamiltonian with Hamiltonian func- 
tions 

H = Hx+aH 2 . 
In the third case, since a — H 2 , we have 

{x\Hx} + H 2 {x\H 2 } = {x\Hx + 

where x l , i = 1,...,6 are coordinates Mi, M 2 , M3, Ti, T 2 , T3. Thus, the system is 
also Hamiltonian with the Hamiltonian function 

Hi 

H = H 1 + ^. 

Regarding integrability of the given five cases, we have the simple Proposition. 
Proposition 5.5. [21] 

(a) A function F is a first integral of equations (15. 4[) if it satisfies 

F={F,Hx} + a{F,H 2 } = 0. 

(b) The Casimir functions Fx and F 2 and functions Hx and H 2 are integrals of 
system (|5.4p for any polynomial a. 

Finally we have: 

Theorem 5.2. |21j System (15. 4p in cases (i)-(iii) is completely integrable in the 
Liouville sense. In cases (iv) and (v), system (15. 4|) is integrable in quadratures. 

5.1. Algebro-geometric integration procedure of the systems. The algebro- 

geometric integration procedure for the first three cases (i)-(iii) (when system (|5.4[) 
is Hamiltonian) is done in |21j . It is based on a construction of the Baker- Akhiezer 
vector- function. 

As usual, we consider the following eigenvalue problem 

4+A(A)W^) = 0, 



dt 

L(\Mt,P) = f&{t,P) 



with a normalization 



* 1 (0,P) + * 2 (0,P) = 1, 
where P = (A, fi) is a point on the spectral curve C. 

Let us denote by oo + and oo _ the two points on the curve C over A = 00, with 
(i = il 2 yj Xq + Zq and \i = —H 2 \/Xfi + Zfi respectively. 



34 



BORISLAV GAJIC 



Proposition 5.6. 21j If the polynomial a(M, T) is a first integral of motion, then 
the vector-function ty(t,P) satisfies the following conditions: 

(a) In the affine part of the curve T, the vector- function ty(t,P) has two time 
independent poles, and each of the components \D (t,P) and ty 2 (t,P) has 
one zero. 

(b) At the points oo + and oo~ , the functions ^f 1 and "J 2 have essential singu- 
larities with the following asymptotics: 



V 1 (t,P) = 



AW X S+ z o (A+a)+7*)t 



+ 0(yj), forP-Kx>- 
■KVXS+ZS (A+a)+^) t(0( i ))) f()r 



* 2 (t,P) 



HV x S+ z o (A+a)+f-)t (O(i)) ) 



for P — > oo 



<W*2+% (*+«)+£> (1 + 0(1)), /orP 
(c) XTie asymptotics have the form 



oo 



f x (t,p) 



^(V^o+^o (A+a) + ^-)t 



1 



# 2 (t,P) = e 



1 



+ 0(1/A 2 ) , P^oo- 



Now we will give explicit formulae for the Baker- Akhiezer function in terms of 
the Jacobi theta-function 9h{z\t) with characteristics [4, i]. 

Let us fix the canonical basis of cycles A and B on T (A ■ B = 1), and let w be 
the holomorphic differential normalized by the conditions 

lu = 2iir, 

I A JB 

A theta-function 6h(z\t) is defined by the relation 



h(z\t) = ^exp 



-r(n + ^) + (z + nr)(n + -) 



Let + and SI be differentials of the second kind with principal parts — | yXo+^o ^ 
and + § v A 2 + dA at oo + and at oo - respectively, normalized by the condition 
that j4-periods are zero. Let us introduce differential 51 = + + f2~. We will denote 
by U the P-period of differential f2, and by c + and c~ the constants: 



fi 



0(1/A), P^P^ 



A 2 + Z 2 A + c- + <3(1/A), P^P- 
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Proposition 5.7. 



-c\ exp 



* 2 (t,P) 



=c 2 exp 



0- 



TTie Baker- Akhiezer functions are given by 

n (A(P + oo+ - Pi - P 2 ) 



i 

—a - 
2 



2 a 



2 V 



2 V* 



■«7) 



6>ii(^(oo+ + oo- - Pi - P 2 ) + <C/) ' 

6> n (-4(P + op- - Pi - P 2 ) + tU) 
9 11 {A{oo+ + oo- - Pi - P 2 ) + «7) ' 



where constants c\ and c 2 are 

= M-AQP - oo + ))gii(^(oo- - Pi))giiM(oo- - P 2 )) 
Cl eu(^(oo-- 00+))^!! (^(P-PO^u^P-Pa)) ' 
giiC^Cg - QQ-))M-4(oo + - P 1 ))9 11 (A(oo+ - P 2 )) 
C2 e 11 (^(oo+-oo-))^ 11 (^(P-P 1 ))e 11 (^(P-P 2 )) ' 

and .4 is the Abel map, and Pi and P 2 are tte poles of the function 4". 

5.2. Classical integration of the systems. In new coordinates 

X-i = olM x + PM 3 X 2 = M 2 X 3 = -@Mi + aM 3 
Yt = oTi + /3T 3 Y 2 = T 2 Y 3 = -0T! + aT 3 

differential equations of motion (|5.4p become 



X 1 = 

x 2 
x 3 

Y = -LV ; V> - X,Y. 
Y 2 

V3 = - \ > - A',V.) - </ v /X : 7 Z,7l' 2 . 

The first integrals are 



XI + Zl (Y 3 + aX 3 ) 

-^XZ + Zi(Y 2 +aX 2 ) 
1 

h 
1 

i 



(5.5) 



Pi 


= XxYt + X 2 Y 2 + 


X 3 Y 3 = ci 


F 2 


= Y? + Y 2 2 + Yi -- 


= c 2 


Hi 


xl + xl + xl 


+ Y l = di 


2I 2 \J Xq + Zq 


H 2 


= X 1 = d 2 . 
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Introducing polar coordinates X 2 = pcosa, X3 — p sin a and using integrals 
(|5.5[) . after simplifying, and denoting p 2 — u, one has 



// :! 



(5.6) u 2 = --J - Bu 2 - Cu - D 

1 < 



2 



where 



B 



h + n 



r -A(X 2 -u 7 2 \( A 2 r 1 d 2 ( Cl -d 2 A) 
C-4(X 0+ Z )(A -c 2 + — 7 === 

D = A(X 2 + Z 2 ){ Cl - d 2 A) 2 

A = di- 



4 



2I 2 \J Xq + Zq 

So, the following proposition is proved: 

Proposition 5.8. The function u(t) is an elliptic function of time. 

Let us remark that u (and consequently p) does not depend on a choice of the 
polynomial a. 

Having u(t) as a known function of time, one can find p(t) as a known function 
of time. In order to reconstruct X 2 and X3, one needs to find a as a function of 
time. 

We have 



1 



p 2 (t) 



x i + z o 



J ( J d 2 +p 2 (t) \ 2 .. 

ci - d 2 di . , ., +ap 2 {t) 

V 2I 2 ^X 2 + Z 2 



The right-hand side of the last equation is a function of time and of the polynomial 
a. When a is a first integral of motion, then the right hand side of the last equation 
is a known function of time. Hence, one can find a by quadratures. In the fourth 

case 

a = XoT, + z r 3 = ^x 2 + z 2 Yl = y/xs + zsdx di tg) 

So, in this case a(t) is a known function of time and one can find a by solving a 
differential equation. Similarly, in the fifth case 

a = Ml + Ml + M| = Xl + Xl +Xl = d 2 2 + p 2 {t) 

is again a known function of time and a differential equation for determining a 
can be solved. Knowing p and a as functions of time, one can easily reconstruct 
X 2 ,X 3 ,Y 1 ,Y 2 andFs- 

Two elliptic curves appeared here. The first one C, has been defined by the 
equation (|5.3p . and it was the curve from which we started. The other one C , 
given by 

3 



(5.7) v 2 = - Bu 2 -Cu- D 
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corresponds to the solution of differential equation (|5.6p . A natural question is how 
these two curves are related. We have the following proposition: 

Proposition 5.9. [21] The elliptic curves C, defined by equation (I5.3[) and C de- 
fined by (|5.7p are isomorphic. 

Using the Sklyanin magic recipe, from Lax representation (|5.ip the separation 
variables for the cases(i)-(iii) are constructed in [21] . Also a sort of separation 
variables for the Hess-Appel'rot case are found. 

6. Motion of rigid body in ideal fluid. Kirchhoff equations 

The mechanical system similar to the motion of a heavy rigid body fixed at a 
point is a motion of a rigid body in an ideal incompressible fluid that is at rest at 
infinity. The equations of motion are derived by Kirchhoff in 1870 (see [37] )• They 
can be written in the form 

(6.1) . d f 

dM 

where Hamiltonian H is homogeneous quadratic function in M and F given by: 
H = l -(AM, M) + (BM, f) + ~<Cf , f). 

Here M is impulsive moment and T is impulsive force. The matrix A is positive- 
definite, the matrices B and C are symmetric. Equations (j6.1[) are Hamiltonian in 
the standard Lie-Poisson structure given by (|3.4[) . Hence for complete intcgrability 
in the Liouville sense one needs one additional independent integral. 

The equations of motion of a heavy rigid body fixed at a point (|3.2I) can be 
written in form (16.11) with H from (13.31) as a Hamiltonian function. 



6.1. Integrable cases. We will list the integrable cases. For a full list and details 
see for example |12) . 

The first nontrivial integrable case of equations (|6.1[) was discovered by Kirchhoff 
in 1870 (see |3Z])- It is defined with conditions: 

• Kirchhoff 's case (1870): 

A = diag(ai,ai,a 3 ), B = diag(6i, b\, 6 3 ), C = diag(ci, c x , c 3 ). 

An additional integral is F4 = M3. It is analogous to the Lagrange case of motion 
of a heavy rigid body fixed at a point. 

• The first Clebsch case (1871): 

A = diag(a, a, a), B = 0, 

The additional integral is: 

F 4 = a(CM,M) -det(C)(C7- 1 f, f). 
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• The second Clebsch case (1871): 

A = diag(ai, a 2 ,a 3 ), 5 = 0, C = diag(ci, c 2 , c 3 ) 
(6.2) c 2 - c 3 c 3 - ci ci - c 2 _ Q 

Oi a 2 a 3 

Conditions (|6.2|) are equivalent to: 

c 2 - c 3 c 3 - Ci ci - c 2 



01(03-03) 02(03-01) 03(01-02) 
where 01,02,03 are pairwise distinct. The additional integral is: 

Fi = 6{M,M) - (Af,f). 

• Stcklov's case (1893): 

^4 = diag(oi, 02, 03), -B = diag(^o 2 a 3 , //a 3 ai, /xaia 2 ), 
C = diag(^ 2 ai(a 2 - a 3 ) 2 , /i 2 a 2 (a 3 - oi) 2 ,^ 2 a 3 (ai - a 2 ) 2 ) 
where /i is a constant. The additional integral is: 
F 4 = Y,( M ! ~ + a* 2 ((o 2 - o 3 ) 2 r 2 + (03 - oi) 2 r 2 + (01 - o 2 ) 2 r 

3 

• Lyapunov's case (1893): 

A = diag(l, 1, 1), B = diag(-2/xdi, -2jud 2 , -2/xds), 
C = diag(M 2 (rf 2 - d 3 ) 2 ;A/ 2 (d 3 - d 1 f,ii 2 {d 1 - d 2 ) 2 ) 
The additional integral is: 

F 4 = d i M j + 2/i(o' 2 d 3 Af 1 r 1 + d 3 diM 2 T 2 + d 1 d 2 M 3 r 3 )+ 
j 

M 2 (rfi(rf 2 - d 3 ) 2 r 2 + d 2 (d 3 - di) 2 r 2 + d 3 (di - d 2 ) 2 r 2 ) 

• Sokolov's case (2001): 

01 = 02 = 1, 03 = 2, 613= a, b 23 = /3, 
C12 = -4a/3, en = Af3 2 , c 22 = 4a 2 , c 33 = -4(a 2 + /3 2 ), 
and additional integral is 

F 4 = (M 3 -ar 1 -/?r 2 ) 2 p + Q 2 , 

where 

P =(a 2 + /3 2 )(M 3 + 2«r! + 2fiT 2 ) 2 + {pMi ~ aM 2 ) 2 
Q = [oAfi + /3Af 2 + (a 2 + /3 2 )r 3 ] (M 3 + 2aTi + 2/?r 2 )+ 
3(/3Mi -aAf 2 )(0Ti -ar 2 ). 

• Chaplygin's first case (1902): 

A = diag(a, a, 2a), £> = 0, C — diag(c, — c, 0), 
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On the symplectic leaf given with (M,L) = 0, the equations admit additional 
integral: 

F A = (Ml - M| + cTlf + 4M^M|. 
• Chaplygin's second case (1897). 

Chaplygin's second case had an invariant relation instead of a fourth integral. 
It is defined in 1897 by Chaplygin (see [13] )■ This system was also considered by 
Kozlov and Onischenko in [IT]. It is defined by: 

A = diag(ai, a 2 ,a 3 ) 





hWa 2 


- a\ T {h 


- MvVb - 


«2 


(6.3) 


bi 3 y/a 3 


- a 2 ± (b 3 


- b 2 )%Ja 2 - 






C\ 3 \/a 2 


— aiT (<=2 


- Ci)y/a S - 


a 2 



0, Ci 2 = 



The invariant relation is: F4 — M\^/a 2 — a\ T M 3 y/a 3 — a 2 = 0. 

Conditions (|6 . 3[) may be regarded as analogy of the Hess-Appel'rot conditions 
in the case of motion of a heavy rigid body fixed at a point. We have shown that 
Hess-Appel'rot case can be considered as a perturbation of the Lagrange top. 

Similarly, the Chaplygin case is a perturbation of the Kirchhoff case. If one 
chooses the basis where a\ — a 2 , the Chaplygin conditions become (see for example 

mna): 

ai = a 2 , 0,13 ^ 0, B = diag(6i, 6i,fe 3 ), C = diag(ci, ci, c 3 ). 

In new coordinates the invariant relation is M3 = 0. 

In the case 5 = 0, Kirchhoff's case can be regarded as a special case of the 
Clebsch case. In [43] Perelomov constructed the Lax representation for the Clebsch 
case as well as higher-dimensional generalizations. Using this Lax representation 
in [35] the Lax representation is constructed for the Chaplygin's second case: 

Theorem 6.1. |22j When B = 0, on the invariant manifold given by the invariant 
relation, the equations of motion of the Chaplygin's second case are equivalent to 
the matrix equation: 

L(X) = [L(X),Q(X)} 
where L(X) = X 2 L 2 + XL\ — Lq, Q(X) = XQi + Qo, and 

L 2 = diag(ci/ai,a/ai,c 3 /ai), Qi = diag(ai, a u a 3 ) 



u = 





M 3 
-M 2 



-M 3 


Mi 



M 2 
-Mi 




U = rr 1 



Qo 





a 3 M 3 + a 13 Mx 
-a x M 2 



-a 3 M 3 - ai 3 Mx 


aiMi + ai 3 M 3 



aiM 2 
-aiMi - ai 3 M 3 
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The spectral curve det(L(X) — ■ 1) = is 

C: n 3 +fi 2 F 3 - A 2 M 2 (c 3 + 2 Cl )+ 

\ln[2F t - (2 Cl + 03)^3] + A^ci( Cl + 2c 3 )- 
\\c\c 3 - A|(2ciFi - d(ci + c 3 )F 3 ) + A 2 ai F 2 = 0, 

where Ai = ^7=. It is singular and has an involution <r : (Ai,/i) — > (— Ai,/x). The 
curve C\ — C/a is a nonsingular genus one curve. 

6.2. Four-dimensional Kirchhoff and Chaplygin cases. In [22] the four-dimen- 
sional generalization of the Kirchhoff and Chaplygin cases is constructed on e(4). 
Let us consider the Hamiltonian equations with Hamiltonian function: 

2H = AijuMijMu +2^2 BijkMijTk + ^ C u T k Ti 

in the standard Lie-Poisson structure on e(4) given by: 

{M, r M,i\ = SikMji + SjiMik - 8 a M jk - 6 jk M a 

{M i:j ,r k } = SikTj - SjkTi 

A four-dimensional Kirchhoff case should have two linear first integrals: Afi 2 and 
Af 34 . It is interesting that under such assumption, the "mixed" term in the Hamil- 
tonian is missing. 

Proposition 6.1. [22] If M\2 and M 34 are the first integrals, then Bijk = 0. 

The proof follows through direct calculations. 
Definition 6.1. The four-dimensional Kirchhoff case is defined by 

2H K =A 1212 M\ 2 + A 1313 (M 2 3 + M 2 4 + Af 2 3 + M| 4 ) + A^M^ 

Ai 234 m 12 M34 + c n (r 2 + r 2 ) + c 33 (rj + r 2 ) 

On e(4) the standard Lie - Poisson structure has two Casimir functions: 

Fx =r 2 + r 2 + r 2 + r 2 , 

f 2 =(M 13 r 4 - M 14 r 3 + Af 34 ri) 2 + (Af 23 ri + M 12 r 3 - A/ 13 r 2 ) 2 + 
(Af 24 ri - M 14 r 2 + M12P4) 2 + (Af 23 r 4 + M 34 r 2 - A/ 24 r 3 ) 2 

consequently, the general symplcctic leaves are 8-dimensional. For complete in- 
tegrability one needs four first integrals in involution. In [22] it is proved that 
except Hamiltonian, the four-dimensional Kirchhoff case has two linear first inte- 
grals F 3 = Afi 2 , F 4 = Af 34 and one additional quadratic first integral: 

F 5 = ai {M 12 M 34 + M 14 M 23 - M 13 M 24 ) 2 

- Cl ((Af 13 r 4 - A/ 14 r 3 + A/ 34 r!) 2 + (Af 23 r 4 + A/ 34 r 2 - Af 24 r 3 ) 2 ) 

- c 3 ((A/ 23 r 1 + A/ 12 r 3 - Afi 3 r 2 ) 2 + (M 24 ri - Afi 4 r 2 + A/ 12 r 4 ) 2 ) 

So, we have 

Theorem 6.2. |22j The four dimensional Kirchhoff case is completely integrable 
in the Liouville sense. 



RIGID BODY DYNAMICS: CLASSICAL AND ALGEBRO-GEOMETRIC INTEGRATION 41 

In the case of the four-dimensional Chaplygin case, one can naturally assume 
that M\2 and M34 are invariant relations. From this assumption, we get: 

Definition 6.2. [22] The four-dimensional Chaplygin case of the Kirchhoff equa- 
tions on e(4) is defined by the Hamiltonian: 

2H Ch =A 1212 Ml 2 + A 13X3 (M? 3 + Ml A + M| 3 + M| 4 ) + A 3434 Ml 4 + 
A 1234 M 12 M 34 + A 1213 M 12 M 13 + A 1214 M 12 M 14 + 
A 1223 M 12 M 23 + A 1224 M 12 M 24 + A 1334 M 13 M 34 + 
A 1434 M 14 M 34 + A 2334 M 23 M 34 + A 2434 M 24 M 34 + 
Bi 2 iMi 2 Ti + B\ 22 M\ 2 Y 2 + Bi 23 Mi 2 T 3 + Bi 24 Mi 2 T 4 + 

-E?341-M34ri + B 342 M 34 T 2 + B 343 M 34 T 3 + £344^/34^4 + 

c u (rf + r 2 2 ) + c 33 (Tl + Tt). 

One can easily check that in this case M\ 2 and M34 are really the invariant 
relations. 
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